home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 19 / CU Amiga Magazine's Super CD-ROM 19 (1998)(EMAP Images)(GB)[!][issue 1998-02].iso / CUCD / Programming / LEDA / prog / graphics / leda.c < prev    next >
C/C++ Source or Header  |  1994-08-05  |  5KB  |  200 lines
  1. #include <LEDA/graph.h>
  2. #include <LEDA/stream.h>
  3. #include <LEDA/array.h>
  4. #include <LEDA/window.h>
  5. #include <math.h>
  6.  
  7.  
  8. typedef GRAPH<vector,int> polyhedron;
  9.  
  10. void rotate(float alpha1,float alpha2, vector& p)
  11. { 
  12.   // rotate 3d-point p about the origin
  13.   // by alpha2 in yz-plane and alpha1 in xy-plane
  14.  
  15.     double R  = hypot(p[1],p[2]);
  16.     double phi = asin(p[1]/R);
  17.   
  18.     if (p[2] < 0) phi = M_PI - phi;
  19.   
  20.     p[1]  = ((R != 0) ? R*sin(phi+alpha2) : 0);
  21.     p[2]  = ((R != 0) ? R*cos(phi+alpha2) : 0);
  22.  
  23.     R   = hypot(p[0],p[2]);
  24.     phi = asin(p[0]/R);
  25.  
  26.     if (p[2] < 0) phi = M_PI - phi;
  27.   
  28.     p[0]  = ((R != 0) ? R*sin(phi+alpha1) : 0);
  29.     p[2]  = ((R != 0) ? R*cos(phi+alpha1) : 0);
  30. }
  31.  
  32.  
  33. inline point project(vector p)   // project p into xy-plane
  34. { return point(p[0],p[1]); }
  35.  
  36.  
  37. void draw_poly(window& W, polyhedron& poly, vector trans)
  38. { edge e;
  39.   forall_edges(e,poly) 
  40.   { point a = project(poly[source(e)] + trans );
  41.     point b = project(poly[target(e)] + trans );
  42.     W.draw_segment(a,b,blue);
  43.    }
  44.  }
  45.  
  46.  
  47. void make_side(polyhedron& poly, node* L, float z)
  48. {
  49.   L[0]= poly.new_node(vector(-70,-20,z));
  50.   L[1]= poly.new_node(vector(-40,-20,z));
  51.   L[2]= poly.new_node(vector(-40,-10,z));
  52.   L[3]= poly.new_node(vector(-60,-10,z));
  53.   L[4]= poly.new_node(vector(-60,+20,z));
  54.   L[5]= poly.new_node(vector(-70,+20,z));
  55.  
  56.   poly.new_edge(L[0],L[5]);
  57.   for(int i = 1; i<=5; i++) poly.new_edge(L[i], L[i-1]);
  58.  
  59.   L[6] = poly.new_node(vector(-30,-20,z));
  60.   L[7] = poly.new_node(vector( -5,-20,z));
  61.   L[8] = poly.new_node(vector( -5,-10,z));
  62.   L[9] = poly.new_node(vector(-20,-10,z));
  63.   L[10]= poly.new_node(vector(-20, -5,z));
  64.   L[11]= poly.new_node(vector( -5, -5,z));
  65.   L[12]= poly.new_node(vector( -5, +5,z));
  66.   L[13]= poly.new_node(vector(-20, +5,z));
  67.   L[14]= poly.new_node(vector(-20,+10,z));
  68.   L[15]= poly.new_node(vector( -5,+10,z));
  69.   L[16]= poly.new_node(vector( -5,+20,z));
  70.   L[17]= poly.new_node(vector(-30,+20,z));
  71.  
  72.   poly.new_edge(L[6],L[17]);
  73.   for(i = 7; i<=17; i++) poly.new_edge(L[i], L[i-1]);
  74.  
  75.   L[18]= poly.new_node(vector( 5,-20,z));
  76.   L[19]= poly.new_node(vector(20,-20,z));
  77.   L[20]= poly.new_node(vector(35,-10,z));
  78.   L[21]= poly.new_node(vector(35,+10,z));
  79.   L[22]= poly.new_node(vector(20,+20,z));
  80.   L[23]= poly.new_node(vector( 5,+20,z));
  81.  
  82.   poly.new_edge(L[18],L[23]);
  83.   for(i = 19; i<=23; i++) poly.new_edge(L[i], L[i-1]);
  84.  
  85.   L[24]= poly.new_node(vector(15,-10,z));
  86.   L[25]= poly.new_node(vector(20,-10,z));
  87.   L[26]= poly.new_node(vector(25, -5,z));
  88.   L[27]= poly.new_node(vector(25, +5,z));
  89.   L[28]= poly.new_node(vector(20,+10,z));
  90.   L[29]= poly.new_node(vector(15,+10,z));
  91.  
  92.   poly.new_edge(L[24],L[29]);
  93.   for(i = 25; i<=29; i++) poly.new_edge(L[i], L[i-1]);
  94.  
  95.   L[30]= poly.new_node(vector(40,-20,z));
  96.   L[31]= poly.new_node(vector(50,-20,z));
  97.   L[32]= poly.new_node(vector(55, -5,z));
  98.   L[33]= poly.new_node(vector(65, -5,z));
  99.   L[34]= poly.new_node(vector(70,-20,z));
  100.   L[35]= poly.new_node(vector(80,-20,z));
  101.   L[36]= poly.new_node(vector(65, 20,z));
  102.   L[37]= poly.new_node(vector(55, 20,z));
  103.  
  104.   poly.new_edge(L[30],L[37]);
  105.   for(i = 31; i<=37; i++) poly.new_edge(L[i], L[i-1]);
  106.  
  107.   L[38]= poly.new_node(vector(55, 0,z));
  108.   L[39]= poly.new_node(vector(65, 0,z));
  109.   L[40]= poly.new_node(vector(60,15,z));
  110.   poly.new_edge(L[38],L[40]);
  111.   poly.new_edge(L[39],L[38]);
  112.   poly.new_edge(L[40],L[39]);
  113.  
  114. }
  115.  
  116.  
  117.  
  118. void make_logo(polyhedron& poly)
  119. { 
  120.   node L[41];
  121.   node R[41];
  122.  
  123.   make_side(poly,L,-5);
  124.   make_side(poly,R,+5);
  125.  
  126.   for(int i = 0; i<41; i++) poly.new_edge(L[i],R[i]);
  127.  
  128. }
  129.  
  130.  
  131.  
  132.   
  133.  
  134. main()
  135. { 
  136.   int w,h;
  137.  
  138.   char* map = Read_Leda_Bitmap("../../incl/LEDA/bitmaps/leda.lbm",w,h);
  139.  
  140.   window W(w,h,window::center,window::center);
  141.   W.set_frame_label(
  142.               "The LEDA Platform for Combinatorial and Geometric Computing");
  143.   W.init(-100,100,-100);
  144.   W.set_show_coordinates(false);
  145.   W.insert_bitmap(w,h,map);
  146.   W.set_node_width(4);
  147.   W.set_mode(xor_mode);
  148.  
  149.   node v;
  150.  
  151.   int speed = 30;
  152.  
  153.   // define a polyhedron in 3d space
  154.  
  155.   polyhedron poly;
  156.  
  157.   make_logo(poly);
  158.  
  159.   //W.draw_point(0,0);
  160.  
  161.   forall_nodes(v,poly) rotate(0.3,0.5,poly[v]);
  162.  
  163.   polyhedron poly0 = poly;
  164.  
  165.   vector dir(2);   // direction and speed of rotation
  166.  
  167.   vector trans0(0,0,0);
  168.   vector trans(0,0,0);
  169.  
  170.   draw_poly(W,poly,trans);
  171.   
  172.   for(;;)
  173.   { dir[0] += random(-10000,10000)/1000000.0;
  174.     dir[1] += random(-10000,10000)/1000000.0;
  175.     dir = dir.norm()*(speed/500.0);
  176.  
  177.     //trans[0] += random(-1000,1000)/5000.0;
  178.     //trans[1] += random(-1000,1000)/5000.0;
  179.     //trans[1] += random(-1000,1000)/5000.0;
  180.  
  181.     if (W.get_button() && W.confirm("quit leda demo ?")) break;
  182.  
  183.     forall_nodes(v,poly) rotate(dir[0],dir[1],poly[v]);
  184.  
  185.     draw_poly(W,poly,trans);   // draw new position
  186.     draw_poly(W,poly0,trans0);  // erase old position  (xor_mode !!)
  187.  
  188.     //trans0 = trans;
  189.  
  190.     // poly0 = poly;
  191.     node w = poly.first_node();
  192.     forall_nodes(v,poly0) 
  193.     { poly0[v] = poly[w];
  194.       w = poly.succ_node(w);
  195.      }
  196.   }
  197.  
  198.  return 0;
  199. }
  200.